STEIN'S METHOD FOR BROWNIAN APPROXIMATIONS 



L. COUTIN AND L. DECREUSEFOND 



Abstract. Motivated by a theorem of Barbour, we revisit some of the classi- 
cal limit theorems in probability from the viewpoint of the Stein method. We 
setup the framework to bound Wasserstein distances between some distribu- 
tions on infinite dimensional spaces. We show that the convergence rate for 
the Poisson approximation of the Brownian motion is as expected proportional 
to A~ 1//2 where A is the intensity of the Poisson process. We also exhibit the 
speed of convergence for the Donsker Theorem and for the linear interpola- 
tion of the Brownian motion. By iterating the procedure, we give Edgeworth 
expansions with precise error bounds. 



1. Introduction 

Among the classics in probability theory, one can cite the approximation in dis- 
tribution of a Brownian motion by a normalized compensated Poisson process of 
intensity going to infinity or the celebrated Donsker theorem which says that a sym- 
metric random walk conveniently normalized also approaches a Brownian motion 
in distribution. Though the topology of the convergence in distribution is known 
to derive from a distance on the space of probability measures, to the best of our 
knowledge, we are aware of only one result precising the speed of convergence in 
one of these two theorems. In [2], Barbour estimated the distance between the 
distribution of a normalized compensated Poisson process of intensity A and the 
distribution of a Brownian motion. The common space on which these two pro- 
cesses are compared is taken as the space of rcll functions, denoted by S([0, 1], R) 
equipped with the distance: 

d (uj,T])= inf (||cjo$ -7^1100 + ||$ -IdHoo), 

$GHom([0,l]) 

where Hom([0, 1]) is the set of increasing homeomorphisms of [0, 1]. It is proved in 
[2] that the speed of convergence is not A -1 / 2 as expected but that there exists a non 
negligible corrective term. This additional term exists because the sample-paths 
of the two processes do not really belong to the same space: Continuous functions 
are a rather special class of rcll functions and sample-paths of Poisson process 
even normalized are never continuous whatever the value of the intensity. Thus 
there is an unavoidable gap between the two kind of trajectories in the considered 
approximation. Actually, the additional term is related to the modulus of continuity 
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of the Brownian motion, i.e. in some sense, it measures the cost to approximate a 
continuous function by a purely discontinuous one. 

We circumvent this problem by considering Poisson and Brownian sample-paths 
as elements of the same space. In fact, the Poisson sample-paths, like the trajec- 
tories of the other processes we are considering in this paper, belong to a much 
smaller space than S([0, 1], R). They all are piecewise differentiable, i.e. of the 
form X)«eN r «(^ — ^n)l[t„. t„ +1 ) where (t n , n > 1) is an increasing sequence of real 
and r„ are differentiable functions. An indicator function is not continuous but it 
has more property than being rcll. In particular, it belongs to 1f$, v for any p > 1 
and any /3 < 1/p (see Appendix A for definition). On the other hand, Brownian tra- 
jectories are (1/2 — e)-H61der continuous so that they belong to Xa, P any p > 1 and 
any < 1/2. Therefore, the natural candidates to support both the distribution of 
piecewise differentiable processes and that of the Brownian motion are the spaces 
T/3,2 for any f3 < 1/2. The original problem is then reduced to the computation 
of the distance between between a given distribution and a Gaussian law on some 
Hilbert space. 

The Stein method is known for a long time to give the speed of convergence of 
many Gaussian approximations (see for instance [5] and references therein). The 
usual approach requires some sort of coupling to derive the pertinent estimates. It 
is only recently that the mixing of Stein approach and Malliavin calculus proved 
its efficiency (see [18] for a thorough analysis of this line of thought): The search 
of ad-hoc couplings in the Stein method is there bypassed by using integration by 
parts formula in the sense of Malliavin calculus. In particular, it has been used 
for approximations of point processes functionals [7, 8, 20]. But to the best of 
our knowledge, up to the notable exception of [2], all these investigations consider 
finite dimensional Gaussian random variables. We here develop the framework for 
a Stein theory on Hilbert spaces. It requires two types of Malliavin gradients : One 
used to characterized the target (Gaussian) measure, one built on the probability 
space of the measure to be compared to the Gaussian measure, used to perform 
the necessary integration by parts. We show that our method is applicable in three 
different settings: Whenever the alea on which the approximate process is built 
upon is either the Poisson space, the Rademacher space or the Wiener space. 

Since [1, 13], it is well known that Stein method can also lead to expansions 
of higher order by pursuing the development. We here generalize this method to 
infinite dimension and give the exact values of the bound coefficients. 

This paper is organized as follows. We first construct the Wiener measure on 
the Besov-Liouville spaces. In Section 3, we transfer this structure to Z 2 (N) for 
the ease of computations. Section 4 is devoted to the development of the abstract 
version of the Stein method for Hilbert valued random variables. In Section 5 
to Section 7, we examplify this general scheme of reasoning successively for the 
Poisson approximation of the Brownian motion, for the linear interpolation of the 
Brownian motion and for the Donsker theorem. For the former two examples, we 
iterate the procedure and derive an Edgeworth expansion with explicit error bounds. 
In Section 8, we show that by a transfer principle, similar results can be obtained 
for other Gaussian processes like the fractional Brownian motion, extending some 
earlier result [9]. Technical proofs are deported to Section 9. 
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2. Gaussian structure on Besov-Liouville spaces 

In order to compare quantitatively the distribution of a piecewise differentiable 
process with that of a Brownian motion, we need to consider a functional space to 
which the sample-paths of both processes belong to. Ordinary Brownian motion 
is known to have sample-paths Holder continuous of any order smaller than 1/2. 
Thus, Theorem [A.l] ensures that its sample-paths belongs to 2ja,oo C Xg,2 for any 
(3 < 1/2. 

Moreover, a simple calculation shows that l[ Oi+00 ) = (e a ). Since for any e > 0, 
2 * s embedded into Hol(e) and since e belongs to the dual of Hol(e), it follows 
that e a belongs to Z*i 2+€ 2 = -^-i/2-e,2- Hence l[ a .+oo) belongs toIi/2- e ,2 f° r anv 
e > 0. This implies that random step functions belong to 2^, 2 for any (3 < 1/2. 
The space of choice may thus be any space 1/3,2 for any (3 < 1/2. The closer to 1/2 
j3 is, the most significant the distance is but the the greater the error bound is. 

To construct the Wiener measure on X/3,2, we start from the Ito-Nisio theorem. 
Let (X n , n > 1) be a sequence of independant centered Gaussian random variables 
of variance 1 defined on a common probability space (O, A, P). Let (e„, n > 1) be 
a complete orthonormal basis of £ 2 ([0, 1]). Then, 

B{t) :=Y,X n ll + {e n ){t) 

n>l 

converges almost-surely for any t G [0, 1]. From [14], we already know that the 
convergence holds uniformly with respect to t and thus that B is continuous. To 
prove that the convergence holds in L 2 (Q; Ip.2), it suffices to show that 

(i) E W^Wh,, = E ||/ Ve»l& = ii/^Hhs < 00. 

n>l n>l 

From [24], we know that I 1 ^ 13 is an Hilbert-Schmidt operator from C 2 into itself if 
and only if 1 - £ > 1/2, i.e. /3 < 1/2. Thus, for ^ < 1/2, the distribution of B 
defines the Wiener measure on 1/3,2- We denote this measure by {ip. Note that (1) 
implies that the embedding from Xi_/3,2 into C 2 is also Hilbert-Schmidt and that 
its Hilbert-Schmidt norm is ||/q+^||hs- F° r latter purposes, we set 

c/j = II/q/IIhs = 2r(1 _ j8) ((1 _ /3)(l/2 - /3)) ' 

according to Theorem A. 3. By the very definition of the scalar product on 1/3,2, 
for r\ G I/3, 2, we have 



[exp(i(r?, w)i 3 2 )] = E P 



exp(i £ / (7 1 1 r' 3 o /-/)»7(s) e n (s) ds X n ) 



= cx p(-^ E (/ 1(/ i r/3 7 o"+ )^( s ) e «w ds ) ) 

= ex P (-i||(7 1 1 _-^o7-/) ?? ||| 2([ o,i ] )) 

= cxp(-i jf (Jq+ " o I\Z /3 )t)(s) r)(s) ds), 
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where f] is the unique element of C 2 such that r\ = Iq+t). Thus, \ip is a Gaussian 
measure on 2^ 2 of covariance operator given by 

This means that 

E w [exp(i{r), w) I/3 , 2 )] = exp(-^(V^j?, ??)i^ 2 ). 

We could thus in principle make all the computations in Xs,2- It turns out 
that we were not able to be explicit in the computations of some traces of some 
involved operators the expressions of which turned to be rather straightforward in 
/ 2 (N). This is why we transfer all the structure to / 2 (N). This is done at no loss 
of generality nor precision since there exists a bijective isometry between Ip 2 and 
l 2 (N). 

3. Gaussian structure on i 2 (N) 

Actually, the canonical isometry is given by the Fourier expansion of the /3-th 
derivative of an element of Zp,2- As is, that would not be explicit enough for the 
computations to come to be tractable. We take benefit from the dual aspect of a 
time indexed point process. On the one hand, as mentioned above, the sample-path 
of a point process is of the form 

!>[*„, i](*) 

n>l 

where (t n , n > 1) is a strictly increasing sequence of reals, all but a finite number 
greater than 1, and thus belongs to 2^ 2 for any /3 < 1/2 as shown above. On the 
other hand, it can be seen as a locally finite point measure defined by 

/e£ 2 ([0,l]) 

n>l 

Said otherwise, we have the following definition. 
Definition 1. For (h,u) € 1.\_p 2 x ^32 

1 h(s)dw s :=< h,I- + \u) > x -_^ j+_ i2 =< I^\h),I-P{uj) > £2 . 

Recall that (e„, n <G N) is a complete orthonormal basis of £ 2 ([0, 1]) and set 
h n = 7^r^(e n ). Then (h n , n e N) is a complete orthonormal basis of X^_ fj2 . 
Consider the map Zp defined by: 

Zp : 2£a — > / 2 (N) 

wi — > (^J h n (s)du(s) = (e n , w)c^([o,i]), n > lj . 

Theorem 3.1. The map Zp is a bijective isometry from 1 p t 2 into l 2 (N). Its inverse 
is given by: 

Z/ : l 2 (N)^Tp, 2 

(a n , n e N) 1 — > ^ a rJo+(e„). 

n>0 
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We thus have the commutative diagram. 



Sp—Z-^oVpoZp 



^,2 i 2 m 

According to the properties of Gaussian measure (see [15]), we have the following 
result. 

Theorem 3.2. Let [ip denote the Wiener measure on 1/3,2- Then Z*p^/3 = mp, 
where mp is the Gaussian measure on i 2 (N) such that for any v <G ^ 2 (N), 

f 1 

/ exp(if.u) dmp(u) = cxp(--Spv.v) 

Ji 2 (N) 2 

with the following notations. 

oo oo 

IMI? 2 (N) = ^2 \x n \ 2 and x.y = ^ x n y n , for all x, y E l 2 (N). 

n—1 n—1 

In view of Theorem [3.1], it is straightforward that the map Sp admits the 
representation: 

Sp : l 2 (N) — > Z 2 (N) 
u = (u n , n e N) i — > (5~](h n , hj) C 2 Uj, n e N). 

By C^(/ 2 (N); X), we denote the space of fc-times Frechet differentiable functions 
from / 2 (N) into an Hilbert space X with bounded derivatives: A function F belongs 
to C b fe (Z 2 (N); X) whenever 

\\ F \\cNi2(N);X)-= SU P ess-sup xei2(N) ||V 0) F(a;)|| x ^2 (N) ®i < oo. 
j=i, ••• , fe 

A simple calculation with characteristic functions shows that a probability measure 
m on Z 2 (N) is equal to mp if and only if for any F e C 2 (/ 2 (N); X), 

[ A !i F{x) dm(x) = where A fi F{x) = x.VF(x) - trace^ o V (2) F(u)). 

Ji 2 (N) 

We are thus led to introduce the Markov process the infinitesimal generator of 
which is A 13 . More precisely, we consider its semi-group, usually called Ornstein- 
Ulhenbeck semi-group as its stationary measure is that of the Brownian motion (in 
fact the image of this measure on Z 2 (N)). 

Definition 2. The Ornstein-Uhlenbeck semi-group on (Z 2 (N), mp) is defined for 
any F G L 2 (/ 2 (N), J", mp; X) by 

PfF(u) = I F{er l u + a/I -e- 2t v) dmp(v), 

Ji 2 (N) 

where the integral is a Bochner integral. It is well known that this semi-group is 
ergodic in the sense that 



P t F(u) t ^° c ) J f dmp, mp - a.s. 
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Moreover, if F £ C£(l 2 (N); X), the function (u i-> P t F(u)) is (k + l)-times differ- 
entiable and we have (see [15]J: 

-(fc+i)t 



(2) V (fc+1) (P t f)(u) = H I (V^F)(e- t u + Vl^lF^v)® vdm fj (v). 



Z 2 (N) 



Lemma 3.1. If F belongs to C^(/ 2 (N); X), then, for any integer k, 

/•oo 

/ ess-sup xe ; 2(N) \\V {k+1) (PtF)(u)\\ l2{N)mk+1) ® x dt < Pk+i\\F\\ C k {l 2 {N) . x) , 

v 

where 

_ (l-2/?)- 1 / 2 (2-2/3)-V 2 k 1 
Pfc " 2r(l-/3) l 2' 2 j ' 

and £> is t/ie Beta function. 

In the sequel, we need to apply A 13 to Hilbert valued functions. We thus introduce 
the notion of partial trace. For any vector space X, Lin(X) is the set of linear 
operator from into itself. For X and Y two Hilbert spaces, the partial trace operator 
along X can be defined as follows: it is the unique linear operator 

trace x : Lin(X <g> Y) — > Lin(Y) 

such that for any R G Lin(Y), for any trace class operator S on X, 

trace x (S <8> R) = tracer (S) R. 

For Hilbert valued functions, we define A 13 as follows. 

Lemma 3.2. The domain of A 13 contains C 2 (7 2 (N); X) and for F in this space, 
we have: 

{A F)(u) = u.(VF)(u) - trac e/ 2 (N) (S>V 2 F(u)), for all x E l 2 (N). 

4. Stein method 
For n* and v* two probability measures on R N , we define a distance by 

p%{v\ pf) = sup / FAv* - ( FAp*. 

\\F\\?<lJ J 

where T is a normed space of test functions (the norm of which is denoted by ||.||<x). 
If % is the set 1-Lipschitz functions on ? 2 (N), then p% corresponds to the optimal 
transportation problem for the cost function c(x, y) = ||a; — 2/|U 2 (rsr) ; x i V <= R- N ( see 
[25]). For technical reasons (as in [21]) mainly due to the infinite dimension, we 
must restrict the space T to smaller subsets. We thus introduce the distances Pj 
for j > 1 as 

Pj{u*, p*) = sup f FAv* - ( Fdp*. 

However, these weaker distances still metrize the space of weak convergence of 
probability measures on l 2 (N). 

Theorem 4.1. Let (v*, n > 1) be a sequence of probability measures on l 2 (N) such 
that some j > 1, 

Pj(" n , M ) > 0. 
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Then, (u*, n > 1) converges weakly to p* in Z 2 (N): 
for any F G C°(l 2 (N), R). 

Proof. As Hilbert spaces admit arbitrarily smooth partition of unity [16], for j > 2, 
one can mimick the proof of [10, page 396] (see also [4]) which corresponds to p\. □ 

Say that p* — mp is our reference measure, that is the measure we want the 
other measures to be compared to. Stein method relies on the characterization of 
mp as the unique measure on / 2 (N) such that any F sufficiently regular, 



X. 



A p F{x) dm(x) = 0. 

/; 2 (N) 

It follows that the Markov process of infinitesimal generator A 13 has mp as sta- 
tionary measure. Since its semi-group (Pf , t > 0) is ergodic, by the well known 
propoperties of Markov processes [11], one can write, for F sufficiently regular, 

/ F{x)dmp{x) -F{x) = [ A^PfF(x) dt, z/*-a.s. 

JP(N) JO 

Thus, 

p 2 (v*,mp) = sup ( [ A P t F(x)dtdu*{x). 

Thanks to the integration by parts induced by Malliavin calculus, we can control the 
right-hand-side integrand and obtain bounds on pi(y* , mp). To be more illustra- 
tive, the Stein method works as follows: construct a process (t M> X(x, t)) constant 
in distribution if its initial condition x is distributed according to mp. Moreover, 
for any initial distribution, the law of X(x, t) tends to mp as t goes to infinity. Stein 
method then consists in going back in time, from infinity to 0, controlling along the 
way the derivative of the changes, yielding a bound on the distance between the 
two initial measures. Other versions (coupling, size-bias, etc) are just other ways to 
construct another process X. In these approaches, for every v* , the couplings are 
ad-hoc whereas Malliavin calculus gives a certain kind of universality as it depends 
only on the underlying alea. Malliavin structures are well established for sequences 
of Bernoulli random variables, Poisson processes, Gaussian processes and several 
other spaces (see [22]). In what follows, we show an example of the machinery for 
each of these three examples. 

The core of the method can be summarized in the following theorem. 

Hypothesis I. For X a Hilbert space, H G Z 2 (N) ® X and a a non-negative real, 
we say that Hyp(A", H, a) holds whenever for any G G C 2 (/ 2 (N); / 2 (N)) 

(3) E„. [Id .G] - E„. [trace(tracejf {H ® H ) o VG)] 

< a \\G\\c%(l 2 (N);l 2 (TSi)®3) \\H\\p(N)®X- 

Theorem 4.2 (Stein method). Assume that Hyp(X, H, a) holds, then, 



(4) P2(v*, mp) < (3 2 trace(traceA-(# ® H) - Sp) 



+ a/3 3 \\H\\% 



(N)®X- 
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Remark 4.1. T/ie iwo terms in the right-hand- side of (4) are of totally different 
nature. The trace term really measures the effect of the approximation scheme 
whereas the second term comes from a sort of curvature of the space on which is 
built the approximate process. As will become evident in the examples below, this 
term is zero when the Malliavin gradient satisfies the chain rule formula and non- 
zero otherwise. 

5. Normal approximation of Poisson processes 
It is well known that for a Poisson process N x of intensity A, the process 

N\t) = -±=(N x (t) - Xt) 

converges in distribution on D to a Brownian motion as A goes to infinity. For any 
f3 < 1/2, we want to precise the rate of convergence for the distance dp, where for 
two measures \x and v supported onI^ 2 , 



dp(n, v) = sup 



11 ciVe 



;R) 



By isometry, this turns out to be equivalent to estimate P2{2pv x , mp) where v x is 
the distribution of N x onl^ 2 - 

Let X[o, l] the space of locally finite measures on [0, 1] equipped with the vague 
topology. We identify a point measure uj 
process 



SneN w i tn the one dimensional 



N 



: t€ [0, 1] .— ► ( du(s) = V l [0 ,t](t n ) 



The measure v\ is the only measure on (x[o, i]i ^(X[o, l])) such that the canonical 
process TV is a Poisson process of intensity A dr. It is customary to define the 
discrete gradient as 

D T F(N) = F(N + e T ) - F(N), for any r e [0, 1], 

where N + e T is the point process TV with an extra atom at time r. We denote by 
D 2j i the set of square integrable functionals F such that E J Q |D T ,F(iV)| 2 dr 
finite. We then have the following relationship: 



is 



(5) 



g(T)(dN(r)-XdT) 



X E 



D T F g{r) dr 



for any g e C 2 ([0, 1]) and any F G D 2 ,i- Moreover, 



It is clear that 



D T (J g(s)(dN(s) - Xds)^j =g(r). 
ZpN x =^-y=J Q h n (s)(dN(s) - Xds), neN 



For the sake of notations, we introduce 



Rx = dfjN x and H x = ^— h v jeN 
so that 3p(Nx + e T ) = 3fjN x + H x (t) = R x + H x (t). 
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Theorem 5.1. We denote by v* x the distribution of2pN\ in Z 2 (N). The measure 
v* x satisfies Hyp(£ 2 ([0, 1]), Hi, l/2\/A). Hence, 



(«) P2^\, rnp) < ^-J= Cp. 



Remark 5.1. It is remarkable that by homogeneity, the partial trace of Hi ® Hi 
is equal to Sp. The only remaining term in Theorem [4-2] comes from the fact that 
the discrete gradient does not satisfy the chain rule. 

One could also remark that the choice of the space in which we embed the Poisson 
and Brownian sample-paths ( i. e. the choice of the value of P) modifies only the 
constant but not the order of convergence, which remains proportional to A -1 / 2 . 

As is clear from the proof, we may expect to obtain a refined approximation if 
we go further in the Taylor expansion of the discrete gradient. By doing this, we 
obtain the so-called Edgeworth expansion. 

In what follows, we make the convention that a sum like Yl°r=i ■ ■ ■ 1& zero - Let 
(a(r, k), r > 0, k > 0) be a family of real numbers such that 

a(0, 0) = 1, a(r, k) = if k - r 2{1, • • • , r} 

and for any r G N, any m G {1, • • • , r} 

1 r 1 

(7) a(r, r + 2m) = r + 2m ^ ( J - + 1 ) ! Qi ( r ~ 3, r-j + 2(m - 1)). 
Consider also the sequence (7^, s > 0) given by the recursion formula: 

s+l ^ c s+3 

( 8 ) 7s = J2 J\ 7s+i-J + ( s + 2)\ /3s+3 " 
Theorem 5.2. Let s be an integer. For F G Cj^ s+1 ^ +1 (Z 2 (N); R), we have 

(9) / F{u)dv* x (u)= { F(u)dmp(u) 

J/ 2 (N) J/ 2 (N) 

+ vf-L) Va(r, r + 2fc) / V (r+2fe) F(u). / Hi{T)® {r+2m) ah dmp(u) 

r^i VvA/ ^ ,A 2 (N) Jo 

+ Rcm(s, F, A) 

w/iere £/ie remainder term can be bounded as 

\Rem(s, F, A)| < 7s A"( s+1 )/ 2 ||F|| c , (s+1)(i2(N);R) . 

Corollary 5.1. Le£ s be an integer. For F G C 2s+1 (7 2 (N); R), we have 



[ F{u)dv* x (u)= [ F(u)dmp(u) 

JZ 2 (N) J/ 2 (N) 

+ V -7= > a(r, r + 2fc)/ F(u)H r+2m( / Hi(r)dT.u) dmp{u) 
7~i V v A/ ^ j; 2 (n) Jo 



+ Rem(s, F, A), 
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where (Hk, k > 0) is the family of Hermite polynomials defined by 

n k (xMx) = (-i) k -^(x), 



with &(x) = cxp(-x 2 /2)/V2^. 

Remark 5.2. In [3, 13], the following polynomials are introduced. For any u£R, 
there exists a family of real polynomials (Q r {-), r > 1) such that 



(10) 



(oo r+2 \ oo 



/n particular, Qo = 1. -For r > 0, it is well known (see [3]) that Q r has degree 3r, 
we denote (q r (k), k = 0, • • • , 3r) its coefficients. For a function F , we introduce 
the formal polynomial Q r {F) defined by substituting t k by J F(u)V.k{u) dmf>{u) in 
the expansion of Q r {t). According to [3, 13], if we choose as functional F{N\) — 
(j)(s) dN\(s) where cf> G C°°([0, 1]; R) and 4>{s) 2 ds = 1, their asymptotic Edge- 
worth expansion coincides with (9) for any order. By identification, it follows that 
a(r, k) — q r (k) for any r and k. Thus, Equation (7) gives a recursive but simple 
method to compute the Edgeworth coefficients up to any order. 

6. Linear interpolation of the Brownian motion 
For m > 1, the linear interpolation B^ of a Brownian motion is denned by 

m— 1 

5^(0) = and dB^i) = m ^ + 1/m) - B\i/m))l [i/mAi+1)/m) (t) dt. 

3=0 

Thus, 3/3 -Bm is given by 

M+Vj/m \ 

m {B\i + 1/m) - B\i/m)) / fc„(i) dt, n > 1 . 

Consider the £ 2 ([0, l])-ortho normal functions 

ef(s) = v / ml [i / Tni(i+1 )/ m )(s), j = 0, •• • , m - 1, s e [0, 1] 

and i 7 ^ = span(e™, j = 0, • • • , m— 1). We denote by p F t the orthogonal projection 
over i 7 ^. Since i?"^ is constructed as a function of a standard Brownian motion, 
we work on the canonical Wiener space (C°([0, 1]; R), Xi,2, m'). The gradient we 
consider, D\ is the derivative of the usual gradient on the Wiener space and the 
integration by parts formula reads as: 



[ u(s)AB\s) =E mt / D\Fu{s)As 
Jo J Uo 



(11) E mt 
for any u G £ 2 ([0, 1]). 

Theorem 6.1. Let be the law of ZpB^ on l 2 (N). The measure satisfies 
Hyp(£ 2 , Hi, 0) where 

Moreover, 
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where 



sup 

-0, 2" 



|Holo(l/2-/3)- 



Remark 6.1. for this example, the constant a of Hyp(X, H 7 a) is zero because 
satisfies the chain rule for derivation. It follows that the Edgeworth expansion 
has a completly different flavor from that established for Poisson approximation. 



Let 



73 



C/3 03 cj 



and 7] 



Theorem 6.2. For any integer s, for any F G C 2,S+2 (7 2 (N); R), we have the 



2 s {s + iy. 

-b 



3/2-/3 

• any in 

following expansion: 

E »l = oTT\ [ (V (2]) F(u), (Si S^'h W 



+ Reni^F, s, m), 

where = trace^ (H^ <g> H^) and Rcm^(F, s, m) can be bounded by 

it \\ F \\cl s + 2 {P(^);-R) 



Rem'(F, s, m) 



< 



m ( S +l)(l/2-/3) 
7. DONSKER THEOREM 

The same approach can be applied to have precise asymptotics for the Donsker 
theorem. Let X = (X n , n e N) be a sequence of independent and identically 
distributed Rademacher random variables, i.e. P(X n = ±1) = 1/2 for any n. For 
any k in N, we set 



, X/.-1, 1, Xk+i • • • ) 

and Xf, = (Xi, • • • , X^-i, — 1, -Xfc+i ■ ■ ■ )■ 
The discrete gradient on this probability space is given by 

D{F{X)= l -{F{X+)-F{X^)). 
Then, the integration by parts formula reads as 

.fceN J L feeN 

for any u = (itk, fceN) which belongs to Z 2 (N). 
Consider 

j [mt] m ^ 

S m(*) = —j= ), X k = z2 —?= X k l[fc/m, 1](*)- 

Since the derivative of with respect to t is a sum of weighted Dirac measures, 
we have 



3 Bl = (£^=x k h n 

k=i Vm 



(fc/m), n > 1). 
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Theorem 7.1. We denote by the distribution of^pB^ on l 2 (N). The measure 
vfn satisfies Hyp(Z 2 (N), H^, 1) where 

Hi = (-L h n (k/m) l [hn] (k), n e N, k e N). 



Furthermore, for any e > 0, t/iere exists mo such that for m > mo, 
p 2 (^ )m/3 )<(l + e)^^m' 3 - 1 /2. 

Remark 7.1. For £/ws example, the Edgeworth expansion would be the sum of the 
previous Edgeworth expansion as neither the trace term nor the constant a are null. 

8. Transfer principle 

For X and Y two Hilbert spaces and a continuous linear map from X to Y. 
Let /i and two probability measures on X and /iy (respectively i/y) their image 
measure with respect to 6. Since is linear and continuous, for F e F e C^(F, R), 
F o 9 belongs to C£(X, R), hence, we have 

sup / F d,uy — / F d^y = sup / f oSdp- / F o 9 d^ 

FeC£(Y,R)J J FeCj (F, R) J J 

< sup / Fdfi- / Fdz/. 

As an application, we can precise the convergence established in [9]. Note that 
in this paper, the key tool was also a matter of Hilbert-Schmidt property of some 
operator. 

The fractional Brownian motion of Hurst index H e [0, 1] may be defined (see 
[6]) by 



B H {t)= f K H (t,s)dB(s), 



where 



The Gauss hyper-geometric function F(a, /3, 7, z) (see [17]) is the analytic continu- 
ation on C x C x C\{-1, -2, . . .} x {z E C, Arg\l - z\ < tt} of the power series 

(a)k(P)k k 

and 

Y(a + k) 

(0)0 = 1 and (a) k = — = a(a + 1) . . . (a + k - 1). 

1 (a) 

Furthermore, according to [23], Kh is a continuous map from L 2 in to 

hence the map Oh = Kh°I^+ can be defined continuously fromX ( g i 2 toIii-(\/2-p), i- 

Since 9#F = B H , we have the following result. 

Theorem 8.1. For any H £ [0, 1], for any 1/2 > e > 0, 



A» 3 

27x 1/2 - £ - 
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9. Proofs 

Proof of Theorem [3.1]. The norms on Ig" 2 and 2 are defined such that 

IMIz+ 2 - IMIz+_ li2 - 

It follows that Ig >2 and 2 are isometrically isomorphic. Since T^_ x 2 is an 
Hilbert space, Parseval equality implies that 



MI I+ 



= V(/i„,cj) 2 _ + , 

-1,2 z ' -^1-/3,2* £ — 1,2 

n>l 



for any CONB (ft n , n > 1) of 2 and Qp is an isometry. The inverse of is 
clearly given by 

9- 1 : I 2 (N)— >Z+ 2 
(a„, n e N) i — >■ ^ 

n>0 

The proof is thus complete. 

Proof of Corollary [3.1]. By a change of variables, it is straightforward that 

di=I B(* 1), 
2 v 2 2 ; 



o Vl - e~ 



2/, 



where £> is the Beta function. Furthermore, since 3/3 is an isometry, according to 
Theorem [3.2], we have 



(N) 



Ml dm^(y) 



< 



L 1 


i 


i 




r(i- 


P) 


i 




r(i- 




i 





I /3 w|| jC 2 ([0il]) d/i /3 (w) 

[\[\t-s)-PdB(s)) 2 dt 
Jo Jo 



E 



1/2- 



(1 -2/3)- 1 / 2 (2-2/3)- 1 / 2 . 



1/2 



The result then follows from (2). 

Proof of Theorem [4.2]. For F e C t 3 , we have 

(12) E„» [F] — E TO(3 [F] 

= - [ [ x.VPfF(x) - trace (s p \7^ pf F(x)) dtdu*(x). 

JP(N)Jo ^ ' 

Applying Hyp(X, H, a) to G = VPfF, we have 



|E„. [F] — E TO/3 [F}\ < 



< 


E„» 


f/ 









+ a||^ll?»(N)®^. 



/ HV^P/Fl^di 
Jo 



< 



/SallVFIUI tracc(trace x (i/ ® ff) - + a&ll# ll?2( N )®xll V (2) i^||oo, 
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according to Lemma [3.1] and to the relation (see [15]) |trace(AoG)| < |trace(A)| ||C||. 
Proof of Theorem [5.1]. For G E C b 2 (Z 2 (N)), we have 



E„ ; [R X .G] = E„. 
According to the Taylor formula 



J (H x (r), D r G( J R A )) ;2(N) Adr 



D T G(R\) = G(R X + H x (t)) - G(R X ) 

= (VG)(R x ).H x (r) + J (l-r)(V 2 G(R x +rH x ),H x (T)®H x (r)) l2{ ^ 2 dr 
Hence, 



E„. [R X .G] = XK* 



J ((VG)(R X ), H x (t) ® # a (t)) ;2(n)82 dr 

(1 - r) (V 2 G(i? A + riJ A ),i7 A (r)® 3 ) /2(N)(9 3 dr dr 



On the one hand, we have 



On the other hand, since V 2 G is bounded, we have 



E„. 



Hence, 



j[ j[ (1 - r) (V 2 G(i? A + ri? A ), ff A (r)® 3 ) i2(N)83 dr dr 



< ^sllV^Hooll^llf. 



(N)®X- 



|E„. [ii A .G] -E„. [trace(VVG)]| < ^||V 2 G|| 00 ||^ 1 ||f 2(N) ^, 

which is Equation (3) with a = \j2\f\. The formula (6) then follows from the 
remarks that trace£2(ifi ® H\) = Sp and that ||77i|| ;2 ( N ) = cp. 



Proof of Theorem [5.2]. We proceed by induction on s. Theorem [5.1] boils 
down to (9) for s = 0. Assume that the results holds for s. According to the Taylor 



16 



L. COUTIN AND L. DECREUSEFOND 



formula, we have 



E V .[F(N X )]= / F(u)dm (u) 

Ji 2 (N) 

s + 2 1 / 1 \ J-l /-oo /■! 



+ E-j(^j 7 i2N y o 7 o V« +1 )P t F(n).(^ 1 (r))®^+ 1 )drdid^( U ) 



J'=2 ' 
X -(s+2)/2 

(s + 2)\ X 

r oo /■! /-I 



/ Iff (1 -^) S+2 V (s+4) P t F(u + 0i/ 1 (r)).(i/ 1 (T))®( s+4 )d0drdtd^( U ) 
J; 2 (n) Jo Jo Jo 

= / F(u)dm^(M) + Ai +^ 2 - 

Ji 2 (N) 



In view of the recursion hypothesis, we have 
3+2 x-0'- 1 )/ 2 s ^'+ 2 / i \ r r 

( 13 ) ^=E^T— E TV E«( r ' r + 2 *) x 

j=2 - y - r=0 \VA/ fe=0 

/ / 1 V^ + 1 + r + 2 ' £ )P t F(w).(if 1 (T))®(-''+ 1 + r + 2fe )dTdtdTO (3 (u) 

J; 2 (N) Jo Jo 

s+2 >_(j_i)/2 /-oo 

+ V / Rem(s + 2-j, V^' +1 )P t F, A)di = Pi+P 2 

~^ J ! Jo 



Now then, according to the very definition of the Ornstein-Uhlenbeck semi-group, 
for any integer I, V^P t F = e~ lt P t V^F. Moreover, since mp is P t invariant, for 
any integer I, 



'Z 2 (N) 



I / V (i) P t F(u).(Pi(r))^ drdidra^ 
o Jo 



u) 



- / / V (Z) F(u).(Pi(r))®Wdrdm^( W ). 

» Ji 2 (N)J0 



We proceed to the change of variables r <— r + j — 1 , j <— j — 1 , m <— k + ( j + 1 ) / 2 
in (13) to obtain 



s+1 / 1 \ r r I / ' 1 \ 
Bi=V -f V tt~ V 7 rMr -j,r-j + 2 (to - 1)) x 



( 2 (N) Jo 



V (r+2m) F(«) . (Pi ( r ))®( r + 2m > dr dm^ (u). 
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Moreover, we have 

s+2 



Rem(s + 1, F, A) = ^ — A - ^ -1 ^ 2 / Rem(s + 2 - j, V (j+1) P t F, A) di + A 2 

s + 2 \-(s+2)/2 roc 

< E J, T.+2-j css-sup„ e;2(N) ^ || V( 2s+7 -^P t P( U )|| dt 



J=2 

\-(s+2)/2 /-l (-oo 

+ (s + 3) , II ^ i?i(r)dr||^ )C s S -sup ue;2(N) ^ ||V( s+4 )F t F( M )||dt 

A -(s+2)/2 

< E 7i 7s+2-j/32 S +7-j||-F 1 || C 2s-3+6 



7 



\-(s+2)/2 f l 

+ 77T3)T 11 i, ff iW d -H^N)^ + 4ll^ll c ^. 



We remark that for s > 0, 



s + 3 < _ max (2s - j + 6) = 2s + 4 = 2((s + 1) + 1). 



Furthermore, 



/W)dr|| 2 2(N) = E (fhn{r)Ar) < ]T [' h 2 n (r) dr = \\c% 



nGN x " u ngN' 

According to Lemma [3.1], it follows that 



Rem(s + 1, F, A) 



< A-^+ 2 )/ 2 f g ± 7hhA^H + 
\i=2 - y - 



3,-1 I || Hi. 



(s + 3)!" s ^; "* llc »" 
In view of (8), the proof is complete. 

Proof of Corollary [5.1]. We need to introduce an auxiliary space, namely 
Z 2 (N) :~ V^ 2 Z 2 (N) equipped with the scalar product 

/ \ T / 1/2 T/-1/2 

\u, «);|(N) = Vp u.Vp v. 

Following [15], the couple (7 2 (N), /^(N)) can be viewed as an abstract Wiener space 
where ^(N) plays the role of the Cameron-Martin space. 

Let denote the gradient on this space: For any Hilbert space X, for any 
v e ^(N), any x in X. any w\, ■ ■ ■ , w n in ? 2 (N), any u e / 2 (N), 

£^/((ioi,u)j2, ••• , (tui,it)j2) <g> a; 

n 

= ^2 d jf(( w i> u )i 2 fi i ■ ■■ . K 1 «)^h®ie lJ(N)«i 

We denote by D2 1 (AT), the completion of the set of such cylindric functions with 
respect to the norm 

11*112,1=/ F 2 (u)dmp(u)+ [ \\D^F(u)\\% { ^ x dm^u). 
J; 2 (N) J/ 2 (N) 
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The n-th power of D 13 and O2, k{X) are denned inductively. Note that if F belongs 
to Cj(Z 2 (N), R), it belongs to the domain of D 2 ,fe(R)- Moreover, it is well known 
(see [19, 22] ) that 



P(N) 



F{u) Hk(u.v) dmp(u) 



for any v € Z|(N). By a density argument, if v € Z 2 (N) and F e B> 2 , fc(R-); then we 



have 

E?^ 

hence the result. 



Z 2 (N) 



F(w) Hk(u.v) dmp(u), 



Proof of Theorem [6.1]. For G sufhcently regular, according to the definition of 
B m and to (11), we have 



(14) 



1,4 [Id.G] =E r , 



m— 1 



" V ~J - p(i-\-l)/m 

E™E + X / m ) _ B{i/m)) / /i n (*) di G n (W m ) 

LnGN i=0 "' i / m 
m-1 r (i+l)/m [ /-(i+l)/m 

/i„(i)dtE mt / D\G n {W m )ds 

Ji/m 

,1 I" /-(i+l)/m 

= / fft,(t)d*.E mt / DlG„(^ m )d S 

Since ft obeys the chain rule formula, 

DtG„ = E V k G n {W m )D\WJ? 



feeN 



m — 1 



= E V fc G„(^ m )(m J] l [l/mt {l+1)/m) (s) / fc fc («) ds) 

fceN ;=o Jl/m 

= VG n {W m ).Hl{s). 

Combining (14) and (15), we get 



E„t [Id.G]=mE t 



m-l Ai+l)/m /-(i+l)/m 

E E E VfeG ™ / h n (t)dt / h k (s)ds 

= E v t m [trace(trace £ 2(ijt i #tj VG)] . 
It follows that i/tj satisfies Hyp(£ 2 , iJ^, 0). To conclude, it remains to estimate 

trace(trace£2 (H^ (g> iftj — S/3). 

In view of the definition ^g, for any integer n, h n is of norm 1 in we thus 

have 

\h n {t) - h s (s)\ <Cp\t- s\ 1/2 - fi , for any t,s£ [0, 1]. 



STEIN'S METHOD FOR BROWNIAN APPROXIMATIONS 



19 



Thus, 

|trace(trace£2(7J^ ® Hi) — Sp)\ 



< 



E 



,-l m-l »(%+!)/ m 

/ h n (s) 2 ds-m^(l K(s)dsf 

Jo -_n Ji/m 



m—1 



■II/. 



(*+l)/m 



-(i+l)/m 

h n (s)(h n (s) - ml 

J it 



(*+l)/m 



A 



Mr) dr) 



ds 



< 



(3/2-/3)mV2-/3 



Ell* 

nGN 



«ll£ 2 - 



The result follows from the definition of cp. 

Proof of Theorem [6.2]. For s = 0, the result boils down to Theorem [6.1]. We 
proceed by induction on s. According to the induction hypothesis, for F sufficiently 
regular, 

/ F(u)dul(u)- f F(u)dmp(u) 

JZ 2 (N) .A 2 (N) 

= [ [°°{VWp t F(u), (Si - Sp)) P{ ^ 2 dtdvl(u) 

Jp(N) JO 

(VW+Vp t F(u), (Si - Sf,)& +2 ) P(N)9V+ , dmp(u) 

+ f f Rem^V^PtF, s, m)dtd^(u). 

Ji 2 (N) Jo 



s 1 

T — 

3=0 J 



2 2 (N) 



Since V^PtF(u) — e Jt P t V^F(u) and since mp is invariant under the action of 
P t , we obtain 



/ F(u)dvl(u) = f F(u)dmp(u) 

+ E 2J+ i C - + 1) , / 2(N) (V(2J+2)F(U) ' {S ™ - S /0® J ' +2 >JW+° dmp(u) 
+ J J™(Rem<< (V™P t F, s, m), (fit, - ^)® 2 ) i 2 (N) « 2 dt d^(«). 

By a change of index in the sum, we obtain the main part of the expansion for the 
rank s + 1. Moreover, 

/>oo 

I Remt(F, s + 1, m)\ < 7 ] / IIV^+^PFIUI^ - S^hw 2 dt 

Jo 

7 t 7 t ll^llc^C+D+Dd^N);^ 



2((s + l) + l) mV2-/3 



hence the result. 
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Proof of Theorem [7.1]. 

E pfcBtG] = ]>>„(fc/™)^G(^i4) 

Let 



.fc=i 



fl-»,(fe) = (-^= h n (k/m), n e N). 

V 777 



According to the Taylor formula, 

= ((VG)(^^), H«,(fc)) I2( N) 
Jo 

+ j[ (1 - r) (V 2 G(^SL - rl4(fc)), 55 ^(fc)) ;2(N)82 dr. 



/;2( N )®2 



dr 



Hence, 



E [a^i.G] = E 



+ E E 

z=±l 



/ (l-r)(V 2 G(a^i + zri/i(fc)),7?i(fc)® 3 ) 

;2( N )®3 ^ r 



According to the definition of V v t , we get 

[3pB* m .G\ -E [trace^aV 2 ^] | < ||V 2 G|UM|? 2(IW(N) , 
which is (3) with a = 1. Furthermore, we have 

|trace(trace;2 (N) (F| l <g> ~ ^/Ol < E( / M s ) 2 ds V] h^k/m) 2 ) 

i=i 7o m fc=i 



+ E 



'»ll£2. 



By the very definition of C^, we have ||/i„||oo < C/3 f° r anv n GN. Thus 



/•l i m m-1 r (k+l)/m 

/ ^ da-- y>(fc/"») 2 < E / "-''^ 2 



_ |/i,(s) 2 -/i J (fc/m) 2 |ds 
fe=o ,/fe / m 

»(fe+l)/m 

<2C^ / /N(fc/m)|da 
fe=o " ,fc / m 

m -i /-(fc+l)/m 

<2( 2 p J2 / Is-fc/ml^ds 
fc=0 " / ' c / m 



2C 



< __^ m /3-l/2. 



3/2- 
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We can choose h„ = V%Iq+ (cos(27rn.)) as a CONB of T\-$,z. According to Par- 
seval formula, 

\\h n \\% = 2 V ( [ h n {s) cos(2irls) ds] . 

1=0 \ Jo / 

Recall that for any (3 > 0, 

r-X 

sup / cos(27rs) ds < oo. 

i£R+ JO 

Since, I 1- ' 9 / can be expressed as a convolution product, there exists c > such 
that 

| / h n (s) cos(2irls) ds\ < c | / s~^ cos(27ds) ds. / cos(27rns) cos(27rZs) ds| 
Jo Jo Jo 

- 1 ! [ s^cos(2^s)ds| l n=l < c f- 1 !^. 
Jo 

Hence, 

]T \\h n \\%<c ^n^-^cm 2 ^ 1 . 

n>m n>m 

On the other hand, 

(m \ 3/2 

v \ fe=o / ieN 



= c , 

lo 



(2-2/3) 3 

Thus, as m goes to infinity, the dominating term is that corresponding to the power 
m P-i/2^ h ence t ne result. 

Appendix A. Besov-Liouville spaces 

This part is devoted to the presentation of the so-called Besov-Liouville spaces. 
A complete exposition can be found in [23]. 

For / e 'C 1 ([0, 1]; dt), (denoted by C 1 for short) the left and right fractional integrals 
of / are defined by : 

Uo + /)(*) = f^j£ /(*)(* - t) a ~ ldt > * > °> 

V?-f)W = f^y J f{t){t-x) a - l dt , x< l, 

where a > and 7° + = I®_ = Id . For any a > 0, p, q > 1, any / G £ p and g e C q 
where + < a + 1, we have : 



f 1 f(s)(IZ + g)(s) ds= 
Jo Jo 



(16) / f(s)(IZ + g)(s)ds= (I?-f)(s)g(s)ds 
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For p e [1, +00], the Besov-Liouville space (£ p ) := I+ p is usually equipped with 
the norm : 

(17) H 7 o+/llz+ p = \\f\\c>- 

Analogously, the Besov-Liouville space (£ p ) :— X~ is usually equipped with 
the norm : 

IIT- a /lk-, P = ll/IU- 

We then have the following continuity results (see [12, 23]) : 

Theorem A.l. i. If < a < 1, 1 < p < then is a bounded operator 

from C p into C q with q = p(l — ap)^ 1 . 
ii. For any < a < 1 and any p > 1, X+ p is continuously embedded in Hol (a — 

1/p) provided that a—1 / p > 0. Holo(^) denotes the space of a Holder- continuous 

functions, null at time 0, equipped with the usual norm. 
Hi. For any < a < f3 < 1, Holo(/3) is compactly embedded in I a ,ca- 
iv. By Iq", respectively I^f, we mean the inverse map of respectively J"_ . 

The relation I^+I^+f = f, respectively I"_lf_f = /, holds whenever 

P > 0, a + (3 > and f e C 1 . 
v. For ap < 1, £/ie spaces X+ and X~ are canonically isomorphic. We will thus 

use the notation X a/p to denote any of this spaces. 

We now recall the definition and properties of Besov-Liouville spaces of negative 
orders. The proofs can be found in [6]. 

Denote by T> + the space of C°° functions defined on [0, 1] and such that 4>^ (0) = 
0, V/c € N. Analogously, set 2?_ the space of C°° functions defined on [0, 1] and such 
that (f>( k \l) — 0, Vfc e N. They are both equipped with the projective topology 
induced by the semi-norms Pk{4>) = S 7 <fcll < A^lloo> Vfc e N. Let T>' + , resp. V_, be 
their strong topological dual. It is straightforward that T> + is stable by and 2?_ 
is stable lf_ , for any j3 e R + . Hence, guided by (16), we can define the fractional 
integral of any distribution (i.e., an element of T>'_ or T>' + ): 

For T € V'_: I^ + T : ^D_4<T, jf_<£ > v >_ >v _, 

For T G /f_T : 0eV+^<T, I§ + <f> >v + ,v + ■ 

We introduce now our Besov-Liouville spaces of negative order as follows. 

Definition 3. For /3 > andr > 1, X~^o r (resp. IZ» r ) is the space of distributions 
T G 2?_ (resp. T € T> + ) such that Iq+T (resp. I^_T ) belongs to C r . The norm of 
an element T in this space is the norm of Iq+T in C r (resp. of lf_T). 

Theorem A. 2. For /3 > and r > 1, the dual space of Xg~ r (resp. Ig~ r ) is 
canonically isometrically isomorphic to I~^ (C r ) (resp. Iq+ (C r ),) where r* — 
r(r-l)- 1 . 

Moreover, for j3 > a > and r > 1, I^_ is continuous from TI ar into I~g~_ a ,.. 

The first part of the next theorem is a deep result which can be found in [24] . We 
complement it by the computation of the Hilbert-Schmidt norm of the canonical 
embedding k 7 from X+ 2 into C 2 . 
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Theorem A. 3. The canonical embedding k 7 from X+ 2 ^ * s Hilbert- Schmidt 
if and only if '7 > 1/2. Moreover, 

1 / 1 \ 1/2 

Proo/. Let (e„, n > 1) be a CONB of £ 2 then (ft„ = Jj+(e n ), n G N) is a CONB 

Of an( i 

ikiil = E = E ii^^iii, = ii^+Hhs 

n n 

and the result follows by straightforward quadratures. □ 
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